Influence of Hydrodynamic Interactions on the Kinetics of Colloidal Particle's 

adsorption 
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The kinetics of irreversible adsorption of spherical particles onto a flat surface is theoretically stud- 
ied. Previous models, in which hydrodynamic interactions were disregarded, predicted a power-law 
behavior t~ 2 ^ 3 for the time dependence of the coverage of the surface near saturation. Experiments, 
however, are in agreement with a power-law behavior of the form t -1 ' 2 . We outline that, when 
hydrodynamic interactions are considered, the assymptotic behavior is found to be compatible with 
the experimental results in a wide region near saturation. 
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The adsorption of colloidal particles or macromolecules 
such as proteins onto adsorbing surfaces is a very com- 
mon phenomenon in many fields of Biology, Chemistry 
and Physics. Deposition of bacteria on teeth or the ad- 
sorption of antibodies on living cells are examples of such 
phenomena which are of great interest in medical sci- 
ences. In many cases, the adsorption is irreversible under 
conditions of practical interest. A model system to study 
the adsorption process is a suspension of latex spheres 
put in contact with a suitable adsorbing flat surface. 
Much work has been done for this system p]-p|,|7]^14[| , both 
on the structural properties of the adsorbed layer and on 
the kinetics of the process. For irreversible adsorption 
of sufficiently light spherical particles onto a flat surface, 
Schaaf et al. M predicted a power-law behavior £ -2 / 3 for 
the time-dependence of the coverage of the surface near 
saturation. This behavior is due to the interplay of the 
diffusion of particles from the bulk and the blocking ef- 
fect caused by the saturation of the surface due to the 
previously adsorbed particles. This result agrees with 
Brownian dynamics simulation of spheres diffusing from 
a bulk solution with a constant diffusion coefficient [[16| . 
Experiments on the adsorption kinetics of small spheri- 
cal particles (proteins or small latex particles) , however, 
show a power law behavior of the form t~ 1 / 2 (T|]2^] when 
gravitational effects on the particles can be ignored in 
front of its pure diffusion. Such a behavior is predicted 
by kinetic models based on the Random Sequential Ad- 
sorption (RSA) filling rules g. RSA, however, ignores 
the physical mechanisms driving the particles to the sur- 
face. Thus, it seems that an important physical ingredi- 
ent has been missed in previous approaches [j. In this 
Letter, we will show that hydrodynamic interactions be- 
tween the particles and the surface substantially modify 
the predicted asymptotic behavior of the system near sat- 
uration. A simple theoretical model taking into account 



the blocking effect, the diffusion of the particles from 
the bulk and the hydrodynamic interactions between the 
particles and the surface, allows us to predict a complex 
asymptotic behavior for the time-dependence of the sur- 
face coverage near saturation. Indeed, near saturation, 
the model predicts a wide first time-domain where the 
time-dependence of the coverage is dominated by the hy- 
drodynamic interactions between the free particles and 
the adsorbing surface. Remarkably, the time-dependence 
in this region is compatible with the experimental find- 
ings of refs. p5| , p5[ . Furthermore, our model also allows 
to derive a second time-domain in which the dynamics 
is dominated by the blocking effect. In this terminal 
regime, the asymptotic time-dependence of the cover- 
age is in agreement with the preditcions of Schaaf et al. 
J7J. Nevertheless, we will see that this terminal regime 
should not be observed for short-range adsorption poten- 
tials, whose interaction range is much smaller than the 
size of the colloidal particles. 

Despite their apparent simplicity, the deposition pro- 
cesses are determined by the interplay of various phenom- 
ena: the Brownian motion of the free particles, the grav- 
itational force, the dynamic interactions mediated by the 
solvent (hydrodynamic interactions) and all other kinds 
of interactions between free particles and the adsorbed 
ones, as well as between the free particles and the wall. Ir- 
reversible adsorption leads to non-equilibrium configura- 
tions, thus, it cannot be studied in the framework of equi- 
librium statistical mecanics. Most of the previous models 
p|,p| JT^JT^ | have neglected the effect of the solvent. They 
consider the particles as moving in dry water [ ]19[ and 
have focused primarily on the geometric aspects, related 
to the excluded surface effects. Recently, however, the 
determinant role played by the hydrodynamic behavior 
of the solvent was pointed out |2Gj2l|] . For instance, the 



theoretically predicted pair distribution function of the 
adsorbed layer shows significant deviations from the ex- 
perimental curves for this function when hydrodynamic 
interactions are ignored [p|, pi)pi| , |2^ , Pq | . The effect of 
the hydrodynamic interaction is to increase the frictional 
force experienced by a particle when it approaches a flat 
surface. Despite its clear implication in the kinetics of 
the adsorption process, the effect of the hydrodynamic 
interaction in the time- dependence of the coverage near 
saturation has never been analyzed before. In the anal- 
ysis proposed here we assume that the free particles dif- 
fusing from the bulk have to cross an entropic barrier, 
due to the presence of the previously adsorbed (bound) 
particles in the layer, before they get trapped by an ad- 
sorbing short range potential between the particles and 
the plane. The density p of free particles in the region 
near the wall is assumed to satisfy a diffusion equation 
of the form §1 
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where 7 = z/2R is the dimcnsionless coordinate in the z- 
direction orthogonal to the wall, R being the radius of the 
particles. The origin of coordinates is taken at the center 
of an adsorbed particle. In this equation, J(-y, t) is the 
flux of free particles in the vicinity of the wall. Thus, it 
is in the region < 7 < 1 that the effect of the excluded 
surface due to the presence of adsorbed particles takes 
place. $(7, 9) is the available area for a free particle to 
move in at a height 7 and at a coverage 9 of the surface 
p4J . 9 is defined as 9 = TrR 2 p s , where the number of 
adsorbed particles per unit of area is denoted by p s . The 
diffusion coefficient D(^,9) is related to the mobility of 
the free particles. Far from the wall, the diffusion coef- 
ficient of spheres in a dilute solution is given simply by 
the Stokes-Einstein formula D = kT/6irr]R, r\ being the 
viscosity of the solvent and is constant. Near the wall, 
however, hydrodynamic interactions modify this behav- 
ior. Lubrication theory J27]] shows that the mobility in 
the direction orthogonal to the wall vanishes linearly with 
the distance between the hydrodynamic surfaces. Thus, 
the diffusion coefficient behaves as 



D(~,,6) = (j + S)D 



(2) 



as 7 — > 0. Do is a constant and 6 = d/2R, where d 
stands for the repulsive range of the adsorption potential 
(fig. 1). Note that the hydrodynamic wall is then shifted 
with respect to the adsorbing wall due to the finite range 
of the adsorption potential considered here. 5 is finite 
but can be made arbitrarily small in our model, in order 
to compare with previous analysis M . If 8 — > no parti- 
cles can be adsorbed in a finite time due to the strength 



of the lubrication forces. The hydrodynamic interaction 
between the free particles and the adsorbed ones is sub- 
dominant for the motion in the direction orthogonal to 
the wall p7] , due to the fact that their surfaces move 
parallel to each other when 7 — > 0. As a consequence, 
the diffusion coefficient is independent of the coverage 9 
in this limit. 

The diffusion equation (|l|) contains all the relevant phe- 
nomena driving the kinetics of particle adsorption. In 
the expression between brackets, the first term stands 
for pure diffusion of the free particles while the second 
accounts for the fact that the available area at a given 
height 7 and at a given coverage 9 is limited by the pres- 
ence of the adsorbed particles. Therefore, if the avail- 
able area is reduced as 7 decreases, this term acts as an 
effective entropic potential tending to decrease the flux 
of free particles. Our model can be applied to a large 
class of systems provided that they meet the following 
requirements; (i) the particles are irreversibly adsorbed 
on the surface and stop moving once trapped; (ii) diffu- 
sion dominates over gravitational effects in the dynamics 
of the free particles. 

Notice that, the use of equation ([j]) implies that we 
describe the transport of the free particles across the 
layer of adsorbed ones without explicit calculation of the 
structural properties of the layer. Brownian motion per- 
mits the particles to explore large regions in space before 
they get adsorbed. Thus, one expects that the overall 
adsorption process is not determined by the local inho- 
mogeneities of the layer of adsorbed particles but by its 
global properties, in the spirit of a mean field approach 
©■ 

The saturation coverage #00, or jamming limit, is 
reached when on the adsorbing surface (7 = 0), the avail- 
able surface function becomes equal to zero. For 9 = 9^, 
the entropic barrier becomes infinite, thus the adsorp- 
tion of new incoming particles is imposible. For spheri- 
cal particles, the entropic potential near saturation can 
be written as 



In $(0, 7 ) -111(600- 9(1 - 7 2 )) ; 



(3) 



The form (9^ — 6>) 3 is the behavior of the available area 
near saturation for irreversible adsorption of particles, 
and has been first derived by Pomeau |2q |. In addi- 
tion, we have explicitly indicated the fact that, at a given 
height, the area excluded by the adsorbed spherical par- 
ticles is reduced by a factor (1 — j 2 ) for spheres P3f . 

In order to describe the kinetics of the adsorption pro- 
cess, we have to find the incoming flux of free particles 
arriving at the adsorbing surface J s (t) = J(j = 0,t). 
Since we are interested only in the kinetics near satura- 
tion, where the adsorption process is very slow due to 



the blocking effect, we can solve eq. (0) neglecting the 
explicit time-dependence of p. We then assume that the 
variations in the density profile and, thus, in the flux, 
adiabatically follow the changes in the coverage through 
$(0, 7) [i||. Therefore, we set dp/dt ~ in the left hand 
side of eq. (Q), implying that J(j, t) is independent of 7 
and equal in fact to J s (t) . We consider here that the den- 
sity of particles in the bulk ps is the control parameter 
and thus express J s in terms of ps, with the boundary 
conditions 



p{l = 1) =Pb, 

p(l = 0) = 



(4) 
(5) 



The first boundary condition assumes that the density of 
bulk particles in the vicinity of the adsorbed layer is ap- 
proximately constant due to the slow adsorption process 
occuring near saturation. We thus consider a particle's 
reservoir located at 7 = 1 with a density p = ps constant. 
The second boundary condition stands for an irreversible 
adsorption: free particles reaching the wall become irre- 
versibly adsorbed and then the density of free particles 
is zero at 7 = 0. We can thus obtain the flux of parti- 
cles reaching the surface in terms of ps by solving the 
differential equation 
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with boundary conditions specified in eqs. (0) and rtq). 
We find the following kinetic equation 
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as 9 — > ^oo. A closed equation for the time-dependece of 
the coverage then follows by multiplying both sides by 
ttR 2 , yielding a generalized Langmuir equation |23| 
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where we have defined the kinetic coefficient K a — Dqtt. 
A crucial point in our analysis is that the leading contri- 
bution to 1(9) near saturation depends on the behavior 
of the integrand for small 7, which allows us to use the 
expression of .0(7) given in eq. (0). Inserting this de- 
pendence in the right hand side of eq. rtq) we obtain that 
the adsorption rate near saturation is proportional to 
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The asymptotics of 1(9) as given by this expression 
strongly depends on the relative magnitude of 5 and 
A9 = (#00 — 9)/9 oc . Clearly, when the coverage ap- 
proaches saturation, in an initial regime the condition 
(6>oo — 9)/9 oc ^> <5 is satisfied since 5 is a constant that can 
be taken as arbitrarily small. In this region, the adsorp- 
tion rate is dominated by the hydrodynamic interactions 
and takes the assymptotic form 



1(9) 
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(11) 



The range of validity of this regime is determined by 
the fact that 1(9) must be positive since eq. (||) de- 
scribes a relaxation process in which the coverage tends 
irreversibily to saturation. Effectively, the condition 
In A9/S 2 -3/2 > indicates that A9/5 2 > exp(3/2) - 1. 
Therefore, the crossover coverage scales as, A9 C ~ 5 2 . In 
this region, the time-dependence of the coverage can be 
obtained by inserting the asymptotic behavior given in 
eq. ( )ll| ) in the right hand side of eq. (O). After integra- 
tion we obtain 
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for A9/S 2 > exp(3/2). Eq. ([l2|) gives an implicit re- 
lation between the time and A9. The scaling of the 
crossover time t c is obtained by inserting the scaling of 
the crossover coverage in this expression, giving t c ~ 
l/<5 4 . Notice the fact that if <5 — s- 0, i c — > 00, indicating 
that this regime must dominate the asymptotic behavior 
of the coverage near saturation. The numerator on the 
right hand side of eq. (Jl3) is a slowly varying function 
of A9. This suggests an iterative procedure to obtain 
the behavior of A9 for times t <C t c . Effectively, one can 
write 
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where, in 



deriving the last expression, a term 
Jhiln(A9/S 2 ) - 1 < ln(l/2S 2 t) has been neglected. Eq. 
(|13| ) predicts a novel behavior for the time-dependence 
of the coverage near saturation dominated by the hy- 
drodynamic interactions between the free particles and 
the wall. Such a behavior differs from that found from 
the RSA model M as well as from that predicted by the 
model incorporating the diffusion of the particles from 
the bulk §. 

The terminal regime t ^> t c or, equivalently, A9/S 2 < 1 
can also be obtained from eq. (|l0|) . The adsorption rate 



in this case is dominated by the blocking effect and obeys 
a different asymptotics of the form 



1(0) ~ 5 (Boo -0) 



-ft) 5 / 2 



(14) 



Notice that the right hand side of this equation vanishes 
as S — ► 0. From cq. (^|) and (|lj) one arrives at the 
A9 ~ £~ 2 / 3 behavior as found by Schaaf et al. J7|. There- 
fore, the model proposed here is also able to reproduce 
Schaaf 's regime M when the adsorption kinetics is dom- 
inated by the blocking effect. However, due to the fact 
that the crossover time between the two regimes scales 
as (5~ 4 and tends to infinity as S — > 0, it suggests that 
this regime is never observed for short range adsorption 
potentials. 

In summary, we have explicitly discussed a model 
where the hydrodynamic interactions are included, in ad- 
dition to other physical mechanisms like diffusion and 
blocking effect, which are relevant to describe the ad- 
sorbing rate. Hydrodynamic interactions play a crucial 
role in the adsorption kinetics and cannot be avoided in 
any experimental work on this process. Previous models 
ignore the physical mechanisms driving the particles to 
the surface (see ref. M and related). Schaaf et al. 0, 
by taking into account the diffusion of the particles from 
the bulk as well as the blocking effect, made a signifi- 
cant step in the description of adsorption kinetics. How- 
ever, the behavior predicted by Schaaf et al. has never 
been observed. Experimental results |lqj2q] suggest a 
behavior near saturation compatible with a power law 
A8 ~ t^ 1 ' 2 . An important conclusion that can be drawn 
from the present work is that the regime predicted by 
Shaaf et al. should not be observed for short range ad- 
sorption potential. The most relevant result of this Letter 
is, however, to have shown that the inclusion of hydro- 
dynamic interactions leads to a behavior near saturation 
compatible with the experimental findings, in view of the 
slow behavior of the logarithmic factor in eq. (|l3|). As 
it has already been pointed out for structural aspects of 
colloidal particles' adsorption onto solid surfaces [|2(j| . pif , 
the results of the present work stress the importance of 
the hydrodynamic interactions also in the kinetics of the 
process. 

P. Wojtaszczyk would like to acknowledge Professor 
J.M. Rubi, G. Gomila and J.M.G. Vilar for fruitful dis- 
cussions and suggestions. This work has been supported 
by the " Ministerio de Educacion y Sciencia" of Spain 
(contract No. SB95-B41700129). 



[1] I. Langmuir, J. Am. Chem. Soc. 38, 221 (1916); 39, 1848 

(1917). 
[2] Z. Adamczyk, M. Zembala, B. Siwek and P. Warszynski, 

J. Colloid Interface Sci. 140, 123 (1990). 
[3] Z. Adamczyk, B. Siwek and M. Zembala, J. Colloid In- 
terface Sci. 151, 351 (1992). 
[4] Z. Adamczyk, B. Siwek, M. Zembala and P. Warszynski 

Langmuir 8, 2605 (1992). 
[5] P. Schaaf and J. Talbot, Phys. Rev. Lett. 62, 175 (1989) 
[6] P. Schaaf and J. Talbot, J. Chem. Phys. 91, 4401 (1989) 
[7] P. Schaaf, A. Johner and J. Talbot, Phys. Rev. Lett. 66 

1603 (1991). 
[8] G. Tarjus and P. Viot, Phys. Rev. Lett. 68, 2354 (1992) 
[9] P. Wojtaszczyk, E.K. Mann, B. Senger, J.C. Voegel and 

P. Schaaf, J. Chem. Phys. 103, 8285 (1995). 
[10] P. Wojtaszczyk, P. Schaaf, B. Senger, M. Zembala and 

J.C. Voegel, J. Chem. Phys. 99, 7198 (1993). 
[11] P. Schaaf, P. Wojtaszczyk, E.K. Mann, B. Senger, J.C. 

Voegel and D. Bedeaux, J. Chem. Phys. 102, 5077 

(1995). 
[12] P. Schaaf, P. Wojtaszczyk, B. Senger, J.C. Voegel and H. 

Reiss, Phys. Rev. E. 51, 4292 (1995). 
[13] B. Senger, P. Schaaf, J.C. Voegel, P. Wojtaszczyk and H. 

Reiss, Proc. Natl. Acad. Sci. USA, 91, 10029 (1994). 
[14] E.K. Mann, P. Wojtaszczyk, B. Senger, J.C. Voegel and 

P. Schaaf, Europhys. Lett. 30, 261 (1995). 
[15] J.J. Ramsden, Phys. Rev. Lett. 71, 295 (1993). 
[16] B. Senger, P. Schaaf, J.C. Voegel, A. Johner, A. Schmitt 

and J. Talbot, J. Chem. Phys. 97, 3813 (1992). 
[17] J. Mackenzie, J. Chem. Phys. 37, 723 (1962). 
[18] R. Jullien and P. Meakin, J. Phys. A. 25, L189 (1992). 
[19] R.P. Feynman, The Feynman Lectures on Physics, 

(Addison- Wesley, Reading 1966). 
[20] F.J. Bafaluy, B. Senger, J.C. Voegel and P. Schaaf, Phys. 

Rev. Lett. 70, 623 (1993). 
[21] I. Pagonabarraga and J.M. Rubi, Phys. Rev. Lett. 73, 

114 (1994). 
[22] I. Pagonabarraga, P. Wojtaszczyk, J.M. Rubi, B. Senger, 

J.C. Voegel and P. Schaaf, J. Chem. Phys. 105, 7815 

(1996). 
[23] P. Wojtaszczyk, J.B. Avalos and M. Rubi, EuroPhysics 

Letters (in press). 
[24] B. Widom, J. Chem. Phys. 44, 3888 (1966). 
[25] C.A. Johnson and A. Lenoff, J. Colloid Interface Sci. 179 

587 (1996). 
[26] Y. Pomeau, J. Phys. A 13, L193 (1980). 
[27] W.B. Russel, D.A. Saville and WR. Schowalter, Col- 
loidal Dispersions, Cambridge University Press (1989). 
[28] P. Wojtaszczyk "Proprietes Statistiques d'Ensembles de 

Particules Adsorbees sur une Surface Solide" Ph.D. The- 
sis, Universite Louis Pasteur, Strasbourg (1995). 



FIGURE CAPTIONS 



Schematic representation of the adsorption process. 
The adsorbing (free) particles diffuse in the bulk and get 
finally adsorbed at a distance d (arbitrarily small) from 



the wall. The origin of the dimensionless coordinate 7 is particles, 
taken in the plane defined by the centers of the adsorbed 
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